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Abstract. This article presents an algebraic topology perspective on the problem of finding a 
complete coverage probability of a domain by a random covering. In particular we obtain a general 
formula for the chance that a collection of finitely many compact connected random sets placed on 
a f-complex X has a union equal to X. The result is derived under certain topological assumptions 
on the shape of the covering sets (the covering ought to be good), but no a priori requirements on 
their distribution. An upper bound for the coverage probability is also obtained as a consequence 
of the concentration inequality. The techniques rely on a formulation of the coverage criteria in 
terms of the Euler characteristic of the nerve complex associated to the random covering. 



1. Introduction 

We consider finite random covering^ of a metric space X, i.e. finite collections of compact 
random sets: A = {Aj}, i = 1, . . . ,n, understood as measurable maps [17_, p. 121] 



where {Q,a,F) is an underlying probability space, and C{X) is the set of compact subsets of X 
topologized by the Hausdorff distance, with the associated Borel algebra o"Borei- The set C{X) U {0} 
is a disjoint union with the point {0}, which plays a role of the empty set. A typical example of an 
(infinite) random covering is a coverage process defined on the Euclidean space. In |17j it is simply 
a sequence of random sets: A = {.^i + ^i, ^2 + ^2, • " " > ^fc + ^fc, • • •}, where {A^} is a fixed family of 
subsets of called grains of the process, and ^ = {^j} random points on M". In applications Ai^s 
are often round balls of a fixed radius, and ^ defines a Poisson process (in which case A is refereed 
to clS 8b Boolean model). 

Problem 1. Given a random covering {Aj}, i = l,...,n of a metric space X, find a complete 
coverage probability: P(X C |A|), where |A| = IJi ^i- 

Reviewing the history of Problem [ij it was first considered by Whitworth, [36] in the basic 
case of a finite collection of independent identically distributed fixed length a arcs {Aj} on a 
unit circumference circle. Much later, Stevens J31j provided a complete answer to the question of 
Whitworth in the form 
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The Steven's result was further improved by Siegel and Hoist [30] where they allowed varying 
lengths for the arcs. In |14] . Flatto obtained an asymptotic expression for coverage as a — )• 0. The 
extension of the circle problem to the 2-sphere 5^ was considered by Moran and Groth [28j , who 
derived an approximation for the probability P(S'^ C |A|), and later Gilbert [16] showed the bounds 

(1 - A)" < P{S^ C |A|) < ^n{n - 1)A(1 - Ar^\ 

where A = (sin is the fraction of the surface of 5^ covered by a spherical a~caps i.e caps of radius 
a. For a G [f ,vr], the explicit expression for P{S'^ C |A|) has been found by Miles [26]. Work in [7] 
provides explicit formulas for the complete coverage probability for a-caps on the m-dimensional 
unit sphere when a G [f ,vr], and upper bounds for a G [0, vr). The literature concerning the 
coverage probability in the asymptotic regimes (where the diameter of grains tends to zero) is vast 
and we only list a small fraction here |15[ [29l [21 [271 [S] ■ Further, the reader may consult the recent 
work in |7] for a more accurate account of the history of Problem [l] 

In this article we focus solely on the case of coverage for 1-dimensional domains X which are 
homeomorphic to finite graphs (specifically, we will rely on the notion of the A-l-complex |18)). 
We denote by dX the set of leaf vertices of X or equivalently the boundary set of a geometric 
realization of X via an embedding in M^. Our random covering A = {Aj} on X, will always be 
finite (z = 1, . . . , n). As already mentioned before, the basic example of a random covering is e- 
balls: {B{^i, e)}, z = 1, . . . , n in an intrinsic metric on X, with centers distributed in an arbitrary 
fashion. We approach the coverage problem by deriving a random complex AA(A) from the usual 
"topological" nerve (c.f. [34]) of realizations of A and considering its Euler characteristic which we 
denote by x(A) =x(AA(A)). 

Let be the set of labeled subcomplexes of the n- vertex full simplex A„. By the labeling we 
understand that every subcomplex s E C„ comes with an indexing of its vertices by numbers from 
1, . . . , n. Elements s, r, k G C„ can be identified, in a non-unique way, with subsets of the power set 



2{iv,"} (^ggg Section 2.1 ). For instance a singleton r = {/} (where I C {1, . . . , n}) labels a face of 
An- By a finite random complex on n vertices we understand an arbitrary discrete probability space 
K = (£„,Pk)- In order to define the random nerve AA(A), one builds a distribution on C„ in a way 
dictated by the usual nerve construction, for instance a probability of a fc-face / = {ii, . . . ,ik+i} 
in €n equals 

PI = P({s e I {/} G s}) = ¥{Ai, n A,, n . . . n A^,^, / 0) . (1.2) 

Identity (1.2) can then be extended from faces to subcomplexes i.e. given s G 
P3=F({rG(J:„ I sCr})=P(V|,}e3{nA^/0})> 

Ps = P(s) = P(V{/}e3{n ^ 0}' V{j}^s{ n A, = 0} )• 

iei jeJ 



(1.3) 



Generally, we will make an underlying assumption of the covering A to be good i.e. to satisfy 
(almost surely) the hypotheses of the Nerve Lemma (c.f. Section [4]). In Proposition 4.7 of Section 



4.2[ we show that a covering A of a 1-complex X is always good if its elements Aj are connected 
and small enough in diameter. The first version of our main theorem is stated below. 
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Theorem 1.1 (Coverage probability for compact connected 1-complexes X, with dX = 0). Let 
A = {Aj}, i = l,...,n be a random good covering of X (with dX = 0). Then, the range of 
X = can be restricted to 

rn = x{X)<x{f^)<ri = m, (1.4) 
and the complete coverage probability equals 

¥{X c\k\)=¥{x{k) = x{X)) 

see:„;xW=x(^) sec„ (1.5) 

N 

where (x) = ^ f (x) Cs,fc (x) , N = m- m, 

k=0 



with ps given in (1.3), and 



Cs,k{x) 



Vk{x) = ^^ hh-'-iN-k, vn{x) = ^1 , (1-6) 

m.<ii<i2<-.<jN-k<rn 

^ ^ {-ly-^^^^-^-^i^^^f^^^^ ^fk >ris), 

ifk< r(s), 

where r^ = r^(s), = r^p(s), ?";^^(s) = — r^p(s) stand for a number of respectively total, 
top and lower: even(odd) dimensional faces of s £ and r(s) denotes a number of all faces. 

If the diameter of Ai is smaller than gC almost surely (where C is a length of the shortest cycle 
in X ), then ps further simplifies as 

Ps = IP(Vij6y(s),i^,{A. n A, / 0} ), (1.7) 

where V{s) is the vertex set of s. 

An extension of the above result to the case of 1-complex X, with no assumptions on dX, is 
provided in Theorem |5.1| of Section [5} One obvious corollary of the above result is the fact that the 
complete coverage probability of any good random covering {Aj} is determined by finitely many 
number^ A computational complexity of deriving F{X C |A|) via the method of Theorem 



1.1 



is not addressed here. However, one may expect that, due to the size of the set (tn (see Section 



2.1), determination of as{x) in (1-5) or the set {s G | xi^) = xi^)} to be double exponentially 
hard in n. Note that the coefficients Usix) independent of the underlying distribution vector 
(ps), therefore once computed for a certain size problem can be reapplied as (ps) changes. The 
vector (ps) can be conveniently estimated numerically (e.g. via the standard maximum likelihood 



estimation, c.f. [23]) but again in the simplest case of Equation (1.7) it is of exponential size: 

n(ri + l) 

2 2 . Therefore, in practical situations one may be interested asymptotic distributions of x(A) 
(as n — )• oo) which would lead to parametric estimators or useful bounds for ¥{X C |A|). Currently 
available results (e.g. in |2H I20j ) concern sparse regimes and they are not applicable, unless we 
allow the diameter of random sets in A = {Aj} to tend to as n — )• oo (as is the case considered 



^which is not obvious when considering vacancy, i.e. volume of X — |A| c.f. [17] 
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e.g. in [l5]). As far as the question of useful bounds for F{X C |A|), as a first step we derive an 
upper bound for the coverage probability, via the concentration inequality [3] in the following 

Theorem 1.2. Let A = {Aj}, i = 1, . . . ,n be a random good coveringof X , then 

where fiQ denotes the expected value of the relative Euler characteristic Xre/(A, Agx) of the random 
pair iM{A),M{A9x)). 

Below we outline our strategy to derive formulas of Theorem The first equality in (1.5) 
is a consequence of basic topological considerations provided in Section |4j the second identity in 
(1.5) follows once we determine the probability distribution of x(A) in terms of (ps)- The required 
steps for this are summarized below for a general case of the Euler characteristic of some random 
complex K = (C„,Pk)- Note that Pk is determined by the vector of probabilities (Ps) representing 
atomic measures Ps = Pk({s}) of elements s in C„ or equivalently a vector (ps) where 

Ps=PK({rG(J:„ I scr}). (1.9) 

Thanks to the obvious identity: ps = Sr-scr^r, vectors (Ps) and {ps) are related by an invertible 
binary matrix B 

(ps)=P-(Ps), (1.10) 

determination of is implicit in the following considerations. 

We will think of random variables on K as real coefficient polynomials in the indeterminates 
e/, where ej = ei{K) will be the indicator functions of faces / of A„. Thus the polynomial ring 
M[e/] = M[{e/}/] is of the main interest. We will denote x, /, Q, P, . . . , the algebraic expressions of 
elements in M[e/], and x(K), /(K), Q(K), P(K) . . .etc., when treated as random variables on K. In 
the paragraph below we review a method to obtain general formulas for moments and distributions 
of discrete random variables of finite range. 

Let [m(Q),m(Q)] C Z, m(Q) < m(Q) denote the integer range of a given Q G M[e/]. Also, let 
N = m(Q) — m(Q) and Wfcj(Q) be coefficients of a certain inverse Vandermonde matrix (see Section 
[3]) given explicitly as follows [12j 

Vkj{Q)= ^ ^, i.]tN-k{j){m{Q),---MQ)), k,j = 0,...,N (1.11) 

where eAr_fc(j)(rn(Q), . . . ,m(Q)) denotes a value of the {N — k — l)-th elementary symmetric poly- 
nomial eAr_fc_i on the range of integers m(Q), . . . ,7n{Q) (with the j-th variable set to zero, where 
variables are indexed from to N). Expansions for powers of Q € M[e/] are given as follows 

Q"" = ^ Cs,k{Q)es, es = JJe/. 
seCn les 

Further, for any s S C„ we set 

Q(s) := Q({e/ = 1\I £ s}). 

Namely, Cl(s) is obtained from Q by substituting e/ = 1 for all I £ s. We also denote by Q(s)(0) 
be the constant coefficient of Q(s) (i.e. the value of Q(s) at 0). Moments and the probability 
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distribution of Q(K) can be now expressed as follows 

^(q'=(K)) = ^(Q(s)(0))'=P3 = E ^s.fcW^'^ 

P(q(K) =m(Q)+i) = Yl Ps=YasAQ)Ps, for je[0,N], 

sG€„;Q(s)(0)=m(q)+j seffn 



(1.12) 



where as,j(Q) = J2k=o'^kj{Q) Cs,kiQ) (see Theorem 

Given a random complex K, a basic example of interest is the number of its d-dimensional faces 

fd= Yl ^i- (1-13) 

{/}eir„;|/|=d+i 

and the Euler characteristic of K: 

X(K) : ((J:„,Pk) 
x(K)(s)=x(s)(0) 



3.3 



in Section 



3.2). 



(1.14) 



By the Euler-Poincare formula (see Equation (2.6), c.f. [18j) we have the following relation between 
dTTTsl) and pni) 



n-l 



(1.15) 



d=0 



Thanks to (1.12) to determine distributions of fdi^) or x(K) in terms of (ps) it suffices to obtain 
expressions for the coefficients Cs±{x)- Let = r^(s), = r^p(s), r^{s) = r^ — ftopi^) denote 
a number of respectively total, top and lower: even(odd) dimensional faces of s G Clearly, 

X(s)(0) =x(s) =r+(s)-r-(s). 

The inclusion-exclusion principle yields the following algebraic identitjj^ for Cs,k{fd), 

rtop{s) 

cs,fc(/d)= Y (-i)i^i-i'i(/d(i)(o))'== 5] (-iy^'-(^^-* 

iee:„;ics i=i 
and the following identity for Cs^kix) (wherj^/c > r{s)) 

Cs,k{x)= E (-1)I^I-I'I(X(1)(0))'== Y (-l)l^l-l'l(r-+(l)-r-(l)) 



i 



(1.16) 



iee:„;ics 



ie£„;lCs 



E E(-i) 



rtop(s)-i-j I top 



(1.17) 



i - j + r 



i=0 j=0 

An extension of these considerations to the case of the relative random Euler characteristic Xreii^^ L) 
is stated in Corollary 3^ of Section [3| 

The article is organized as follows: In Section [2] we further discuss the general setup of random 
complexes and their associated invariants, mainly x(A) and associated random variables. In Section 
[3] we derive identities in (1.12), in particular relevant formulas in the case of the random relative 



■^where Cs,k{fd) vanish, unless s is a subcomplex of the complex which is a union of d-dimensional faces, with 
some d-faces present 

^Cs.fc(x) vanishes for k < r. 
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Euler characteristic, in Corollary 3.4 Further, in Section|4]we prove several basic topological results 
showing that the Euler characteristic of the nerve of a covering suffices to determine complete 
coverage of a 1-complex. In Proposition |4.7[ we also provide a sufficient condition for the to be 
good (in terms of a diameter of its elements) and determine when its nerve is equivalent to an 
associated Vietoris-Rips complex. We collect relevant facts in Section [5] and provide the proof of 
Theorem 1 1 . 1 1 there . The upper bound for F{X C |A|) of Theorem |1.2| is shown in Section[6j Section 
[7] concludes the paper with a comparison of our formulas to the Steven's formula As a test 

case we exhibit specific computations in the case of three random arcs. 

Beyond 1-complexes the techniques of algebraic topology provide coverage criteria for higher 
dimensional objects. For instance if an underlying space X is an m-dimensional manifold a necessary 
and sufficient condition for coverage is nonvanishing of the m-th(top) Betti number of the nerve 
AA(A). We aim to develop these ideas in subsequent papers. 
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problems in sensor networks [SI ITUl [5] , and inspiring conversations during the First National Forum 
of Young Topologists at Tulane University in 2009. Both authors acknowledge the support of NSF 
DMS #1043009. The first author was also partially supported by DARPA YFA N66001-11-1-4132. 



2. Random complexes and their topological invariants. 

2.1. Random complexes. We refer the reader to [l5j for background on algebraic topology. 
Consider A„ to be a full simplex on n-vertices indexed by 1, . . . , n (geometrically A„ is the convex 
hull of n points given by the standard basis vectors in M""*"^). Recall that a d-dimensional face in 
can indexed by the collection of its d + 1 vertices: / = {^l,^2^ • • • where 1 < ii < 12 < 

. . . < id+i < n. Denote the set of all faces of A„ by f{n) and particularly d-dimensional faces by 
fd{n), i.e. 

fin) = {I\IC 2ii'-'"}}, fain) = {I\IC 2^''-'^\ \I\=d+ 1} (2.1) 

Next, we want to consider the set of all labeled sub- complexes <tn of A„ union a special point {0} 
with a role of the empty set. By a labeled sub-complex we understand a subcomplex of A„ with the 
labeling given by vertices of A„ . A natural set to consider for enumerating labeled subcomplexes 
is the power set 2/W of /(n) union {0}, which we denote by Here and thereafter, we slightly 
abuse the notation and use s, r, k for elements of both (tn and 

Clearly, there is a surjective correspondence 11 : 1— )■ C„ which to a given subset s G assigns 
a subcomplex n(s) in C„ given by the union of elements {/} of s treated as faces of A^. 11 is 
clearly not bijective, however, with certain choices we may easily build right inverses. In particular, 
we will be interested in two cases, which we refer to as the antichain and chain representations 
^ : <tn ^ ^n, ~ ■ <tn ^ ^n- The antichain representative s G of s G contains only it's top 
dimensional faces, i.e. for any {/}, { J} G s: I J and J ^ / (in terms of the partial order of 
given by c). The chain representative s of s G is obtained from the antichain representative 
by adding all remaining subfaces of s. Clearly, s is the poset associated to the subcomplex s (c.f. 
[33]). We denote the images of C„ in *p„ under ~ and by and respectively. We also have the 
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associated projections 

U-.^n^ ^n, n : ^ (2.2) 
Note that the cardinahty of C^i and therefore C„ and is given by the Dedekind number M(n), 
c.f. [22]. For any s G (tn, we call the elements of s, top faces of s. 

• For the rest of the paper we make the convention to identify €n with C„ i.e. we think about 
subcomplexes as represented by chains (unless indicated otherwise). 

Recall from Section [T| that a finite random complex on n vertices we understand a discrete prob- 
ability space K = (C„,IPk)- It is easy to see that Pk satisfies the following equivalent conditioni 
(for /' C /) 

(A) Pk(/| (I'r) =IPk({sGC„|{I}Gs} I {rGCnlU'j^r}) =0, 

(B) Pk(/' I /) =IP({s G C„|{/'} G s} I {r G £n I {/} G r}) = 1. 

In short (A) says that if a subface /' of / has not occurred then / cannot occur as well, equivalently 
(B) says that /' occurs whenever / has occurred. We say that K is supported on a subcomplex, 
k G if and only if for any {/} ^ k we have Pk({-^}) = 0- Given random complexes K and L on 
n-vertices, the joint probability space (K, L) := (C„ x £„,Pk,l) is a random pair if and only if L is 
almost surely a subcomplex of K, i.e. the following condition holds 

(C) for every (s, r) G €„, x C„ such that r s we have Pk,l( s, r ) = 0. 

For a given K (or (K, L)) it will be convenient to consider Bernoulli random variables which are 
indicator functions of faces in K, i.e. for / G f{n) we define 

ej:€n^{0,l}, ejis) = r' ^' (2.3) 

10, otherwise. 

For any s G Cm we set Cg = H/es indicator function of the subcomplex s, clearly Cg takes 



value 1 on r if and only if s C r. Recalling, vectors (ps) and (Pg) from (1.10), we have the following 
expressions in terms of e/'s: 



(2-4) 

P3=p(( n n (1-^^)) 

lefin),les J6/(n),J^s 



It is clear that one may now construct an inverse of B from (1.10) by expanding the polynomial 

Ps= n «^ n (1-^-^) (2.5) 

/e/(n),/es Je/(n),j^s 



in monomials Cg and taking the expected value (see Section 3.2). 

Summarizing, any random complex K is fully determined by indicator functions {e/} of faces and 
their joint distribution. In the next section ejs will serve as formal indeterminates for functions 
defining topological random variables on K, such as the Euler characteristic. The main use of 



conditions (A) ((B)) and (C) will be to define (in Section 2.3) an appropriate polynomial ring 
where the random topological invariants, such as x(K), live. 



''because events {s G £n 1 {1} £ s} and {r G £„ | {/'} ^ r} are disjoint 
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Remark 2.1. In general, one could consider more flexible model of a finite random complex with 
(*P„,]P) as the underlying probability space. It can be thought of as a distribution on open faces 
of A„ (i.e. interiors of faces) with an exception of the zero dimension (the vertices). In this model 
it is possible, for instance, for an edge to occur without its vertices (i.e. (A) can be violated). We 
have the following obvious result 

Lemma 2.2. Suppose that (^„,IP) satisfies (A) or (B), (with ^„ in place of <tn and¥ in place of 
Pk)- Then P is supported on C„ = C„ C CLnd (Cri,PL ) defines a finite random complex. 

Proof. Let k G it suffices to see that if k 7^ k then P(k) = 0. Pick {/'} € k — k, such that there 
is {/} G k satisfying /' C I. Then k G {s G q3„ | {/} G s} n {r e q}„ | {/'} ^ r} from (A) 

P(k)<P({sGq3„,|{/}Gs}n{rG^?„|{7'}^r}) =0. □ 

Remark 2.3. We may easily generalize the definition of the random complex K to the case n = 00, 
and thus removing dependence in n in the definition. This is done by considering all labeled 
subcomplexes £00 of the infinite simplex Aqq = IJn^"' ^^"^ regarding a random complex K as a 
probability space (£001 Pk)- Such random complex is finite provided the support of K is contained 
in A„ for some n. 

2.2. Topological invariants in the random setting. Recall that, thanks to the Poincare-Euler 
formula [18j, the Euler characteristic of a general n-complex K is given by 

n 

x(i^) = 5](-l)Mim,(i^), (2.6) 

j=0 

where dim.j{K) denotes the dimension of the real coefficient jth chain group Cj{K;W), and equals 
(in the absolute case) to the number of j-dimensional faces fj{K) of K. We will also need a relative 
version of x- Given a pair (K, L) where L is a subcomplex of K we have 

n 

Xrei{K,L) = ^{-iydimj{K,L), (2.7) 
i=o 

where dimj{K, L) denotes the rank the jth relative chain group Cj{K, L;M) = Cj{K;M) /Cj{L;M.) 
(c.f. dS]). Note that we have 

dim,(K,L)=/,(K)-/,(L). (2.8) 

Invariants x{^) Xrei{K, L) can be expressed in terms of the Betti numbers {/3fc(i^)}, {/3fc(K, L)} 
of the chain complexes C*(K) and C*{K,L), (c.f. [18j). Specifically, 

n n 

X{K) = Y,i-^yPj{K), Xrei{K,L) = ^{-iyf]j{K,L). (2.9) 

j=0 j=0 

2.3. Random polynomials. Given a random complex K let us treat the indicator functions of 
faces {e/} (or in a case of a random pair {ej,wj}) as formal indeterminates and consider a poly- 
nomial ring in ej (without loss of generality we work over M): 

]R[e/] := M[e{i}, . . . , e{„}, £{1^2}) • • • 1 • • • 1 ^{i,...,n}]^ 

(or M[e/,'u;j] in case of random pairs). Observe that any random variable X on K is given as such 
polynomial, namely X = ^see:„ ^i^) where Pg is given in (2.5). 
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Based on (2.6) we may express the random Euler characteristic x from (1.15) as the following 



polynomial in R[e/]: 



(2.10) 

/e/(n) 



Lemma 2.4. Given a random complex K = (C„,Pk) o,nd its collection of the indicator functions 
{e/}. Consider Q,Q' G M[e/] as two representatives of the same coset in M[e/]/X, where I is an 
ideal generated by the following relations 

{ejei = ej \ for alll C J}, (2.11) 

(in particular: e| = ej.) Then Q(K) = Q'(K) almost surely. 

Proof. It suffices to show that P(ej ej = ej) = 1 for any /, J where I C J. We have 

P(ej ei = 0)= P(ej = 0, e/ = 1) + P(ej = 1, e/ = 0) + P(ej = 0, e/ = 0). 
Thanks to (A) : P(ej = 1, e/ = 0) = 0, thus 

P(ej e/ = 0) = P(ej = 0, e/ = 1) + P(ej = 0, e/ = 0) = P(ej = 0), 
and P(ej e/ = 1) = 1 - P(ej e/ = 0) = 1 - P(ej = 0) = P(ej = 1). □ 

We will further denote the quotient ring M.[ei]/Z by Mx[e/]. Clearly, Mj[e7] has an additive basis 
of monomials indexed by the chain representatives: s G C^: 

es = llei. (2.12) 

les 

In the case of pairs (K, L) we have a pair of sets of face indicator functions {ei,wj} corresponding to 
K and L respectively. Then, it is relevant to consider a polynomial ring M.[ei,wj] modulo relations 
in (2.11) and additionally (thanks to property (C)): 

{wj wj = wj \ for all / C J}, 
{wi = wjej, I for ail J C l\. 

The resulting quotient ring will be denoted by Mx[ei ,wj], and the analogous statement as Lemma 



2.4 is true for random variables expressed as representatives in Mx[ei ,wj]. An important for us 



example of a polynomial in Mx[e/,ti;j] is the relative Euler characteristic 

XreKK,L) : (e:„ X e:„,PK) ^z, 

Xreii^, s') = Xreii^, s'), if s C s, i.e. the relative Euler characteristic of (s, s') (2-14) 
= 0, if s' s. 

Note, that thanks to (C), the set of pairs (s, s') such that s' ^ s is of measure zero in (K, L) thus 
the value of Xreii^, L) on such pairs is irrelevant. Thanks to (2.8), the polynomial expression for 

Xrel 

(K, L) is given as follows 

Xrei= J]] (-1)1^1^^(6/ -u;/), ^2.15) 
/e/(n) 
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Remark 2.5. One may restrict considerations to the invariant ring Mx[e/]®" or M[e/]®" where (3„ 
is a symmetric group. Indeed, 6„ acts on the vector space V{ei) spanned by {e/ | / G fin)} by 
permuting the faces /, i.e. 

7r(e/) = e/ o vr = e^(/), vr G 6„. 

Recall that M^e/]®" = {Q G M*[e/] | Q = Q o yr}. Noether's Theorem, cf. ^ p. 27], implies that 
e/]®" is generated as a polynomial algebra by {eg | s G £n}, where 



°7r(s) 



and e^(s) = e^f^^j^j^^^j^y^ = e{^(7j)^^(/2)^...^^(7^)}, for s = {h, h, ■ ■ ■ , h}- Note that x G IR*^-/]®'' 
has an expansion 



Hi} 



n 



Hm + 



'{1:2,3} 



+ (-ir 



n 



n 



3{1,...,„}- 



3. Moments and distributions of the random Euler characteristic. 

We begin with basic review of the method of moments for the finite range discrete random 
variable X, and provide a specific formulation based on the recent work in [12]. Alternatively, one 
could use fractional moments (see e.g. [H p. 17]), however they are not easier to consider in the 
setting of the Euler characteristic. 

3.1. Method of moments. First, we need basic information on the Vandermonde matrix V (cf. 
[25]). Given a sequence of values x = {xq,xi, . . . ,xn}, V is an (A^ + 1) x {N + 1) explicitly given 
as follows 

/I Xq ■■■ x^ \ 



V = V(x) 



1 Xi 



V 1 XN ■■■ J 



V is invertible provided the Xj's are distinct (cf. [25]). A closed form of has derived in [T2] in 
terms of the elementary symmetric polynomials. Denote by ei(j)(x) the ith-elementary symmetric 
polynomial in variables: xq, • • • ,Xj, - ■ ■ , xn for j = 0, • • • , N, where xj means that Xj is omitted, 
specifically 

1 if i = 

xi^xi^.-.xi. ifi>0. (3-1) 



e.(i)(x) 



I, l<h<l2<...<h<Ndkf^j 



By [121 p. 647], we have 



V(x) = (vkiix)), where Vkii^) = (-1) 



.N+k ejv,fc(0(x) 

-nN I 



(3.2) 



for z = 0, . . . , A^, A: = 0, . . . , A^. In the case x is an integer interval [m, . . . , m], m, m G Z, m < m of 
size N = rn — m we obtain 



Vki{x) = Vki{m,m) 



(-1) 
A^! 



i+k 



N 



tN~k{i){rn, ...,m), 



(3.3) 
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Lemma 3.1. Let X be a discrete random variable of a finite range x = {xq,xi, . . . ,xn}, and let 
Hk = IE(x'^) denote the k-th moment of X. Given the vector /i = {fiQ, . . . ,fiN) we can recover the 
distribution of X explicitly as follows 



Pi = P(X = Xi) 



N 

E 

k=0 



Vki l^ki 



0,...,iV, 



(3.4) 



where v^i = i'fci(x) are the Vandermonde coefficients. 
Proof. By definition we have a linear system of equations 

N 

fJ'k = ^XiPi, for A; = 0,l,...,iV. 

i=0 

In the matrix form this system of reads: pV = /x where p = (po, . . . ,Pn), and /x = {^q, . . . , ^n)- 
Since all Xj's are distinct det(V) = ni^j(^« ~ ^j) 0- Thus V is invertible and we have the unique 
solution p = fj,V^^. Identity (3.4) is a direct consequence of (3.2). □ 

Our goal for the next subsection is to provide expressions for distributions of polynomial random 
variables in Mx[e/]. 

3.2. Distributions of random polynomials. Since the differences between Mx[e/] and Mx[ei, wj] 



are mostly notational, we will work with the former. Recall from Section 2.3 that any representative 

(3.5) 



polynomial in M[e/] is a linear combination of monomials Ck from (2.12) 

Q = ^ CkCk, Ck G M, 

where the constant coefficient cq = cq is indexed by the empty set. Note that if Q G Mx[e/] then, 
thanks to the relations in Mj [e/] , we may always pick expansions of Q in terms of the antichain or 
chain representatives i.e. 



Q 



or 



Q 



see;„ 



E 



6c; 



(3.6) 



where in the second expansion we just applied our convention from Section 2.1 to identify elements 
of <tn with their chain representatives. We refer to |3.6p eft) as the antichain representative and 
3.6Vight) as the chain representative of Q in Mx[e/]. From the perspective of Lemma 2.4 it is 
irrelevant which expansion of Q we choose. Below, we outline a strategy to determine coefficients 
Ck of (3.5) via the inclusion-exclusion principle. 

Recall, the general form of the inclusion- exclusion principle, |24j : Given a finite set F and 



functions f,g:2 



we have 



satisfying 



S:SCF 



S:SCF 

Recall the following notation: given Q G M[e/] and s G define 

Q(s) := Q({e/ = 1 | / G s}). 



(3.7) 
(3.8) 
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i.e. Q(s) is a polynomial obtained from Q by substituting ej = 1 for all {/} G s, and Q(s)(0) its 
constant coefficient. 



Lemma 3.2. Consider any representative Q G Mx[e/] in a general form (3.5). For any k G *p„ the 
coefficient Ck of Q in the expansion (|3.5[) is ^iuen as follows 



(3.9) 



Ck(Q)= E (-l)W-l^lQ(r)(0). 



/n i/ze case Q is represented by the chain expansion ( right) (2>.&), for any s G s 7^ {0} we have 

Cs(Q)= 5; (-l)N-W(Q(r)(0)-co), (3.10) 

where cq = cq = Q(0) is t/ie constant term ofQ. 

Proof. In the inclusion-exclusion principle set F = k, then any subset S G F is just a subset of 
faces r of k, i.e. r G and r C k. Directly from (|3.5p for any r C k, we have 



Q(r)(0) = c,, 

r'Cr 



thus setting g{r) = Q(r)(0) and /(r) = Cr, Equation (3.9) follows from (3.8). To obtain (3.10) 



consider the polynomial Q = Q — cq. If r C k anqjr 7^ r, then Q(r)(0) = 0. Therefore, for s G £n 



Equation (3.9) yields 



cM= E (-l)l^l-l^lQ(r)(0)= E (-l)'^'-'^'Q(r)(0). 

re?P„,rCs rGCn.rCs 



Because Cs(Q) = Cs(Q) for s 7^ 0, the identity in (3.10) follows. 



□ 



For a polynomial random variable Q G M[e/] in a general form ( |3.5[ ), define constants 

wM='Yc~, = min{cs,0}, m(Q) = ^ ^j", c+ = max{cs, 0}. (3.11) 



Denote the coefficients of the general expansion (3.5) of the chain representative of the fc-th power 
(Q)'' by Cs,fe(Q), i.e. 



Q'' = E Cs,fc(Q)es 



(3.12) 



We summarize efforts of this section by stating the following result which is a direct consequence 
of Lemma 3.1 and Lemma 3.2 and justifies the identity in (1.12) of Section [T] 



"i.e. r is not a set of chains 
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Theorem 3.3. Given Q as a chain representative in Mx[ej], suppose that the set of realizations of 
Q is in the integer interval [rn,fn], then the distribution of Q and its moments are given as follows 

/ifc = E(q'=) = (Q(s)(0))'Ps = Yl Cs,fc(Q)Ps, 

P(Q = m + j)= Y Ps = asj(Q)Ps, iG[0,iV], N = rn-m 

se£„;Q{s)(0)=m+j seC™ 
N 

for as,j(Q) = ^Vfcj(Q)cs,fc(Q), 



(3.13) 



k=0 



where Vkj{Q) were defined in (3.2). Further, cq = Q(0) and Co^k = Cq, and for s ^ 0: 



Cs,fc(Q) 



ree:„;rcs 



|s|-|r| 



(Q(r)(0)-co)^ 



(3.14) 



Proof. Since Cg are Bernoulli random variables 

//fe = E(Q^') = Y Cs,fc(Q)E(es) = Y Cs,fc(Q)Ps 



thus (3.13) is an immediate consequence of (3.4). Formula (3.14) follows from (3.10) applied to 

□ 

3.3. Formulas for x(K), /d(K) and Xreii^: L). In this section we aim to provide a slightly more 
tractable formulas for the coefficients Cs,fc( • ) and the integer ranges [m( • ), m( • )] for the polynomials 
X(K), /rf(K) and 

Xreii^i L), where K is a given random complex on n vertices. These formulas were 
partially listed as (1.16), (1.17) in Section[T] Thanks to Theorem 3.3 it will provide us with a more 
precise characterization of distributions for these polynomials. 

We begin with the case of /d(K). Clearly, the range of fa is contained in between 



in{fd) = 0, and m{fd) 



n 

^d+1, 

For a subcomplex s £ and its corresponding antichain s, recall the following notation 



(3.15) 



^top 


= ^topi^) 


= {numer of 


even dimensional faces in s}, 




r 

top 


= rtopi^) 


= {numer of 


odd dimensional faces in s}, 




^low 


= ^L(s) 


= {numer of 


even dimensional faces in s - 




r 

low 


= ^Fowi^) 


= {numer of 


odd dimensional faces in s — 




'^top 


= rtopis) 


top top 






I'low 


= riowis) 


= s — S , 


r = r(s) = rtop + riow = |s 





(3.16) 



Proof of Formula (1.16). Applying (3.14) directly to fd = /d(K) we obtain the first identity in 
(1.16). For the second equation in (1.16), let 1 E *p„ be the set of all d-faces. Since fd = X]{/}ei 
for any k C 1, Equation (3.9) implies 

c.((/.)')= (-l)l'^l-l'^l(/d(r)(0))^ = J](-1)W-'(IJ (3.17) 

re'P„;rCk i=l ^ ^ 
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Considering f ^ as an element of Mx[e/] and choosing a chain representative for we conclude 
that its coefficients Cs,a;(/^) vanish unless the corresponding antichain s consists of purely d-faces. 
In the latter case we obtain from ( |3.17 ) 



for k = s, 



which implies the identity in (1.16) via the notation of (3.16). 



□ 



Next, we turn to the random polynomial x(K). The range of is contained in [zzi(x); 
where 



E 

r;0<2r+l<n 



n 

2r + 1 



(3.18) 



and 'm{x) = 

r;0<2r<n 

If K is supported on some subcomplex k G smaller than the full n~simplex, the above range 
can be narrowed to 



m.{xm = - /2r+i(k), m{xm= •^^^(k) 

0<2r+l<dim(k) 0<2r<dim{k) 



Proof of Formula (1.17). Applying (3.14) to Q = x directly one obtains the following formula for 
the coefficients Cs,fc(x) 

Cs,k{x)= (-1)I^I-I'I(X(1)(0))'== Y (-l)l^l-l'lK(l)-r-(l))'. 



lGt„;lCs 



lGt„;lCs 



which shows the first part of (1.17). To obtain the second part of (1.17) we can follow the same 
reasoning as for fd above. However, we choose to present a different, slightly more involved, 
argument for the purpose of cross verification. It will use the inclusion-exclusion principle but in 
a different way than in Lemma 3.2, namely via the well known multinomial formula. Recall that 
given indeterminates xi, . . . , Xm we have ^13) 



{Xi + X2 + . . . + Xr^ 



E 



Q;=(«l,a2v,am), 
|Q!|=fc 



^ai a2 c 



(3.19) 



a\\a2^....ar, 



f, CKj > 0, \cx.\ = and cx form all possible partitions of k. Let a have 



where (^) 

coordinates indexed by /(n) (i.e. the faces of A„). A direct application of (3.19) to (2.10) yields 

_E C) n ((-1)"'-'% 



a={ai), 
|q:|=A: 



E ((-i)='"-""-'"-')U>*" 



where we denoted 



Observe that for any a and a'. 



a={ai), 
\a\=k 



s{a) = {le fin) I ai > 0}. 



(3.20) 



-s{a) ^s{a') ) 



m 



(3.21) 
(3.22) 
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if and only if the corresponding antichains are the same i.e. s(a) = s(a'). Fix a chain representative 
of some complex s G £n and let s be the corresponding antichain. Clearly, s C s, consider partitions 
a of k which are in the form a = /3 + 7 where ^ = (/?/), satisfies: f3j > for I £s and /3/ = for 
/ G s — s, and 7 = (77) satisfies: 7/ > for / G s — s and 77 = for I £s. The following claim 
immediately follows 

Claim: Given s G and any partition a oi k indexed by /(n), we have n(s(a)) = s if and only 
if a has the above decomposition: /0 + 7, (see (2.2) for the definition of H). 

Therefore, the Cs^kix) coefficient of the chain representative of (x)^ is a sum of coefficients of 6^(^01) 
for all a in the form +7. Applying notation (3.16) we may express it as 



ix)'' = J2 Cs,fc(x)es, where 

(_l)EliT(l^d-i)ft+E;LT(l^.l-i)7. 

(^,7)=(/3i,...,/3^t„p,7i,...,7r;^^), 
|^| + |7|=fc,ft>0,7j>0 

0, 



(3.23) 



k 



if A; > r, 



otherwise, 



where we indexed the faces of s in s by {li}, i = 1, . . . ,rtop and faces of s — s in s by {Jj}, 
j = 1, . . . ,riow To set up the inclusion-exclusion principle, note that the sum for Cs^kix) is a part 
of the larger sum (where we allow /3j > 0): 



E 



(_l)E;iT(l^d-i)ft+E;LT(l^.l-i)7. 



3lvi/3rtopi7li---)7r;o„), 

|^| + |7|=fc,ft>0,7j>0 



k 



^(_i)(i/d-i) + ^(_r 



(k.l-i) 



i=l 



We stratify the above sum with respect to number of /3j's strictly greater than zero, and set up the 
inclusion-exclusion as follows. Let F = {1, . . . , rtop} and define for any S C F, functions /, g (in 
^7l\) , ([3^) as 

fiS) = 



E 



(/8,7)=({ft},{7i}),l/8| + l7l=fc,7,>0 
/3i>0, if i e 5, ft=0 if j ^ 5. 



ies i=i 
Observe that X^Jif^ = '^low ~ ^/ot«' which yields 

g(S) = (\S+\-\S-\+r+^- 



^low j ' 



where |5'*'|(|S' |) denotes number of even(odd) dimensional faces of s indexed by 5. By (3.8) we 
obtain 

k 



f{F)= Y {-iy''"^-\^\{\s+\-\s-\+' 



S:SCF 
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Since there are even dimensional faces and r^^^ odd dimensional faces in s, for a fixed z G [0, r. 



top 



' topi 



and j G [0,?^^^^] there are exactly ('~*°'') (^*°*') subsets S C F satisfying i = jS""^!, j = \S \. Thus the 



second part of (1.17) now follows from f{F) = Cs,fc(x). □ 
As the last case of interest we consider is the relative Euler characteristic Xre/(K, L) of a random 



pair (K, L). Denoting the characteristic functions of K by {e/} and of L by {wj}, ( |2.7[ ) and (2.8) 
imply the following polynomial expression 

n-l 

Analogously as in the absolute case the distribution of (K, L) is determined by 

Ps,r = IP(es = l,Wr = l)= P(esWr = 1). 

The maximal constants for the range of Xreii^, L) are 

EkiXrel) = mix) - rn{x), and m{xrel) = rn{x) - mix)- 
For convenience we state the following corollary of Theorem |3.3t 



(3.24) 

(3.25) 
(3.26) 



Corollary 3.4 (Distribution oixreli^, L)). Given a random pair (K, L), the distribution ofxreli^, L) 
on [m{Xrei),m{Xrei)] is given as follows, for j £ [0, N], N = m{xrei) - miXrei) 

P(X^e,(K, L) =m(XreO +J) = ^ ^s.r J (Xre/) Ps,r , (3.27) 

{s,r)ee;„xe;„ 

N 

as,T,j{Xrel) = '^[vkjiXrel) Cs,T,k{Xrel)) , 
k=0 



where (using the notation of (3.16)j 

E((XreKK,L))'=) = Cs,r,fc(Xrez)Ps,r, 



(3.28) 



Cs,r,fc(Xrei) — * 



J2 i^_iytopis)+rtopir)-i-j-i'-j' f^topi^)\ f^top(^)\ f^topi^n /'■topW 

ie[0,r+^(s)]jG[0,r-^{s)], V W V W V * / V J 

i'G[0,r+p(r)]j'e[0,r-„p(r)] 



•((i-i)+(*'-iO+(r+„{s)-r-„,(B))+{r+„{r)-r-„(r))) , for k > r, 



0, for k < r. 



The proof is as fully analogous the previous arguments and is omitted. Note that the expression 
for Cs^r,k{Xrei) in (3.28) simplifies to (1.17) whenever L = 0. 
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4. Topological considerations. 



Given a deterministic covering of a finite simplicial 
complex X, i.e. a collection of compact connected 
subsets A = {Ai}, we can define its nerve, M{A) as a 
finite complex where vertices are just elements Ai of 
the covering and a fc-face I = {ii, . . . , ik+i} belongs 
to J\f{A), if and only if Ai^ n Ai^ n . . . n Ai^^^ / 
(c.f. [34J). The following result, due to Borsuk [6], 
is of fundamental importance in algebraic topology 

Lemma 4.1 (The Nerve Lemma [Ql). Let A = {Ai} 
he a covering of X and N{A) the associated nerve. 
If all intersections Ai^ n Ai^ n . . . n Ai^^^ , for k > 
are contractihle, then N{A) has a homotopy type of 
the subspace \A\ = IJi of X. 




Figure [4] An example of A-l-complex with marked 
realization of a good cover. 



If ^ = {Ai} satisfies the assumption of this lemma 
then we call it a good covering (of X). Recall that 
a subset of X is contractible if it can be deformed 

continuously to a point [TH]. In the rest of this section we collect elementary facts from algebraic 
topology and show how the Euler characteristic of AA(A) provides a criteria for a good covering 
A = {Ai}, to completely cover a connected 1-complex X. 



4.1. Coverage via the nerve. We assume throughout that X is a A-complex (c.f. |18l p. 
103]) which is a slight generalization of the notion of the simplicial complex, and in particular 
is more appropriate for quotient constructions. Formally, X is equipped with finitely many vertices 
{vi, . . . Vk} C X and a collection of maps (edges) Sij : [0, 1] X such that Sij{0) = Vi, Sij{l) = vj 



and s 



ij\{0,l) 



is injective, and each point in X — {vi, . . . , Vk} is the image of exactly one such s 



1(0,1) • 



A set C X is open in X if and only if each s^-^(S') is open in [0, 1]. The boundary dX of X is 
the set of leaf vertices of X, when we treat X as a graph (or equivalently the topological boundary 
of a geometric realization of X via an embedding in M'^). 



Proposition 4.2. Suppose X is a finite connected 1-complex, and let {A^} he a good covering of 

(4.1) 



X. Let \ A\ = \JiAi, denote U = \A\ and V = \A\''. Then, 

f3i{X) > /3i([/), 
and 

x{x) < xiu). 



Moreover, if the inequality in (4.1) is strict then (4.2) is also strict. 



Proof. Consider the Mayer- Vietoris sequence applied to U and V: 

^ Fi(c/ n y) 4 Hi{u) e Hi{v) Hi{x) Hq{u n y) ^ Hq{u) e Hq{v) Ho{x) 

Since U DV = OA is just finitely many points, in real coefficients we have 



(4.2) 



0. 
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From xiX) = l-/3i(X), x{U) = Po{U)-pi{U), xiV) = Po{V)-pi{V). Since di is injective 

we have /3i([/) + < which imphes -/3i(X) + /3i(C/) < 0. This proves (O). 



Now to prove ( |4.2[ ) we have two cases to consider: /3o{U) > 1 and /3q(U) = 1. First assume 
Po{U) > 1. We argue by contradiction. That is, suppose xiU) < xi^)- Then PoiU) — f3i{U) < 
(3o{X) - f3i{X) so that /3o(f/) < (3i{U) + 1 - /3i(X). But /3i{A) - (3i{X) < by the previous 
lemma. Therefore we obtain I3q{U) < 1 contrary to our assumption. Now assume Po{U) = 1. Then 
X{U) = 1 - l3i{U) and xiX) = 1 - which yields xiU) - x{X) = -pi{U) + > 0. Thus 

X{U) > x{X). □ 

Suppose A = {Ai\ \s a good covering of X, and |^| = Uj^«) ^y the Nerve Lemma an obvious 
necessary condition for X C is 

x{X) = x{\A\)=x{M{A)). (4.3) 
Clearly, the condition is not sufficient in particular when dX ^ 0, but we have the following 
Corollary 4.3. Let X satisfy dX = then (OJ) implies X c\A\. 



Proof. Notice that generally X (even with dX ^ 0) is homotopy equivalent to a bouquet of circles. 
If 1^1^^ 7^ in X , then (since \A\'^ is open) we pick p £ \A\'^ which is not a vertex of X. Then p 
is in the interior of one of the edges which we denote by e. We may homotopy X away from the 
interior of e to a bouquet of r circles S = \/^ in such a way that p is away from the wedge point 



(just collapse along the edges different from e). From Proposition 4.2 

r-l 

l3i{\A\) < S' V {S' - M)) < /3i(5) = /3i(X). 

Thus /3i(|^|) < /3i{X) and therefore xiX) < xd^l), which implies the claim. □ 



When dX ^ 0, (4.3) is insufficient however we may adjust it by using the relative version 



Xrei{X,dX) of the Euler characteristic in ( |2.7| ) for the pair {X,dX) reduces to 

Xrel{X,dX)=x{X)-#{dX}, 

where i^{dX} is the number of boundary points of X. Because X is a special case of the cellular 
complex, from \18\ p. 102] we have 

Xrel{X,dX)=x{X/dX), 

thus we can work with the quotient X' = X/dX which is a A-l-complex ( |18| p. 103]) with 
dX' = 0. Let q : X ^ X' he the quotient projection, then the covering A oi X projects to the 
covering A' of X' . It is not true that A' is automatically a good covering of X' , one may easily find 
examples where this is the case. However, the following fact is available (easy to see proof omitted) 

Lemma 4.4. Given A = {Ai} is a good covering of X, let for every i the intersection Ai n dX he 
either empty or a point (in other words Aqx = {Ai n dX} is a good covering of dX). Then the 
quotient covering A' of X' is also good. 

Consequently, we say that j4 is a good covering of the pair (X, dX) provided A is good for X and 
Aqx is good for dX. Then by he above lemma A' is good fo X' and Corollary 4.3 says that A' 
covers X', if and only if x(l^'l) = x{X'). To obtain the condition analogous to (4.3), for the 
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general case, we need to express x(l^'l) in terms of the data of A and X. We have the following 



generalization of Corollary 4.3 



Lemma 4.5. Given a good covering A = {Ai} of {X,dX) let \ A\ = {J^Ai. Then X C |^|, if and 
only if 

Xrel{Af{A),Af{Aax)) = Xrel{X, SX) (4.4) 

or equivalently 

X{\A\) = x{X) - ndX} + #{\A\ n dX}. (4.5) 

Proof. Observe that X C \A'\ to X C \A\. Indeed, since |^| is closed if X — 1^41 7^ then we 
may choose a point in x G X — such that x dX, since the projection g is a homeomorphism 



on X — dX, we conclude that q{x) ^ X' — \A'\. Next, Equation (4.4) follows immediately from 



Corollary |4.3[ the fact that A and Aqx are good and the identities 

X{\A'\)=Xreli\A\,\A\ndX), x{X')=Xrel{X,dX). 



Now, thanks to (2.7) we compute 

Xrel{X,dX)=x{X)-#{dX}, 

Xrei{\A\,\A\ n dX) = x{\A\) - #{\A\ n dx}, 



which yields (4.5). □ 



Remark 4.6. Equivalently, the coverage condition for {X, dX) can be obtained by looking at the 
enlarged covering ^, which is A together with the boundary vertices: dX = {xi, . . . , x^jg^}}. 
Then A is good if satisfies the conditions of Lemma 



4.4 



x(\A\) = xm u dx) = x{\A\) + x{dX) - x{\A\ n dx) 
= x{\A\) + #{dx}-#{\A\r^dx}, 



which together with (4.3) leads us to (4.4) 



4.2. Coverage of X by e-balls. Vietoris Rips complex. A special case of interest (see e.g. 
|35t llOj ) is when a connected 1-complex X ought to be covered by e-size neighborhoods, where e 
is allowed to be sufficiently small. In such cases the topology of the nerve of the covering simplifies 
and one may work with Vietoris-Rips complex [19], as we investigate in the following paragraphs. 
First, equip X with an intrinsic distance 

dx(x,y) = min length(7). (4-6) 

7(0)=x, 7(l)=s/. 

i.e. dx{x,y) is the length of the shortest path between x an y, which in practice is just a smallest 
sum of edge-lengths (and their pieces) connecting x and y, (the lengths come from some choice of 
geometric realization of X in M^). We will denote the intrinsic diameter of a subset Y C X hy 
diam(y). Recall that given a simplicial complex K its Vietoris-Rips complex Tl{K) , [19] is defined 
to be a maximal simplicial complex (with respect to inclusion) which has the same 1-skeleton as 
K. Practically, this means that T^iK) is obtained by filling every A;-clique in the graph K^^^ with 
a (A; — l)-dimensional face, e.g. 3-cycles are filled with 2-simplices in TZ{K), etc. 
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Figure 4.2: Possible shapes of closed e-balls in X with the intrinsic distance dx- 

We will consider a finite covering A = {^i, . . . , An] of (X, dx) by closed e-balls, possible shapes 
of such balls for e sufficiently small is depicted on Figure |4.2[ 

Let us denote by 'R-iA) the Vietoris-Rips complex of the nerve of the cover. We record the following 

Proposition 4.7. Suppose C is the length of the shortest cycle in the quotient complex X' = X/dX, 

{i) if e < \C then the covering A by e-balls in {X,dx) is a good cover, 
(ii) if e < then the nerve ^f{A) of A equals Tl{A). 

Proof. For (i) we must show that every /c-fold intersection Ai^ n Ai^ n . . . H Aj^. has a homotopy type 
of a point. Because diam(74i) < C, Ai is a connected tree and therfore contractible, which shows 
the claim for k = 1. For k = 2, first suppose that a nonempty intersection Ai n Aj is disconnected 
i.e. dim(ffo(^j ^ ^j)) ^ 1 (where • ) denotes the reduced homology groups c.f. [IE])- Since Ai 
and Aj are connected, the reduced Mayer- Vietoris sequence for Ai D Aj then simplifies to 

Hi{A^ U Aj) %{Ai n Aj) HoiAi) HoiAj) ^ {0}, 

We obtain Hi{Ai [J Aj) ^ HoiAi n Aj) ^ ]R*= for some k > 1, which implies that Ai U Aj contains 
a nontrivial cycle. This however contradicts the fact that diam(Aj U Aj) < 4e < C. Thus k has 
to vanish and Ai n Aj must be connected contain no cycle and therefore contractible. Now, for an 
induction step with respect to k, it suffices to apply the previous step to A' = Ai^ n Ai^ n . . . n Ai^ 
and A" = Ai^^^ . 

Before proving (ii), recall the 1-dimensional version of Helly's Theorem (c.f. [llj) implies 
that given a finite collection of intervals {Ci, C2, . . . , C„} on R, if the intersection of each pair 
is nonempty, i.e. Ci n Cj 7^ 0, for every 1 < i, j < n, then HILi ^« 7^ ^■ 

First consider the case of 3-fold intersections, i.e. supposing that Aj r\ Aj. ^ 0, 1 < fe 7^ j < 3, 
we aim to show that Ai n ^2 H ^3 / 0. Observe that ^ = ^1 U ^2 U A3 is connected and by 
the argument of (i) it must be a connected tree, i.e. contains no cycles. Let pi^2-, P2,3, Pi, 3 be 
distinct points in V such that pi_j £ Ai D Aj. Note that for each pair: pij, ps^t there exists a 
path in V connecting these points. We now consider two cases: (1) one of these paths, we denote 
by I, contains all three points pij, then the collection {Cj}, Ci = I D Ai, i = 1,2,3 satisfies the 
assumptions of Helly's Theorem which implies the claim. (2) none of the paths between paris of 
Pij 's contain the third point. Consider two shortest paths: li between pi^2 and p2,3, and I2 between 
pi^2 and p2,3 then /i_2 = h I2 is a segment between pi^2 and some vertex of v £ V. The vertex v 
has to be in one of Aj's, w.l.o.g. suppose v £ A2 (as other cases are analogous). Then if v is also 
in Ai or A3 we can take pi^2 or p2,3 equal to v and use (1). If f ^ and v ^ A3 then we observe 
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that either Ai or is disconnected which is not the case. This concludes the proof of (ii) for the 
3-fold case, the general case can be obtained by induction. □ 

5. Complete coverage probability. 
In this section we derive consequences results of Sections |4.1| - |4.2| in the random setting. 

5.1. Random coverings and the random nerve. Suppose A = {Aj} is a random covering of a 
metric space X, we define the nerve M{A) of A by defining a probability measure Pa on <tn via the 
process elucidated in Section [T] in ( |1.2[ ) and (1.3). Observe that given a subspace Y C X we obtain 
an induced random covering Ay from A: 

Ay = {Ainy,A2ny,...,A,ny} 

The definition of Pa extends to pairs (A/'(A), A/'(Ay)) in an obvious way. In particular given (s, r) E 
X we set 

Ps,r = P({(k, 1) G e:„ X e:„ I s c k, r c i}) 

= ip(V{/}6s{n A« ^ 0} '^{j}eA n Ai n >^ / 0} )• ^^-^^ 

i€i jeJ 

Clearly, A^(A) is a random complex, and (A/'(A), A/'(Ay)) is a random pair. We say a finite random 
covering {Aj}j=i^,,.^„ of X is good if and only if it is a good covering on X almost surely. Further, 
we say a random covering A = {Aj} of a pair {X, dX) is good provided it is a good covering of X 
and IKqx is a good covering of dX. |A| will denote the random set Uj ^j. 

5.2. Proof of the extended version of Theorem Let Xre/(A,Agx) be the relative Euler 
characteristic of the pair (A/'(A), A/'(Agx))- We may now state Theorem |1.1| for a general 1-complex 
X. 

Theorem 5.1 (Coverage probability of a 1-complex X with dX ^ 0). Let A = {Aj}, z = 1, . . . , n 
he a random good covering of the pair {X^dX). Then, the range of Xrei{f^i f^dx) can he restricted 
to 

m = Xrei{X,dX) < XreKA, Aax) < u = Tn, (5.2) 
and the complete coverage prohahility equals 

P(X C \k\)=¥{Xrel{A,l^dx)=Xrel{X,dX)), 

asAXrel)Ps,r, ^^'^^ 

(s,r)ec„xe;„ 

where as,r(Xre/) = fls,r,o(XreO are defined in (3.27) of Corollary 3.4. and ps,T in (5.1). 



Proof. Under the given assumptions. Lemma |4.5| implies 

P(X C |A|) = P(XreKA, Aax) = XreliX,dX)). 



(5.4) 



At this point the formula (3.27) of Corollary |3 .4| can be applied to the random pair (AA(A), AA(Agx)) 
to give an exact expression for P(Xre/(A, Agx) = Xreii^, dX)^ In this particular case the range 
of x(A, Agx) is given by (5.2), where the lower bound follows from Proposition 4.2, and the upper 
bound corresponds to the case when elements of the covering A are pairwise disjoint and contained 
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in X — dX, i.e. AA(A) is just n distinct points. The formula for ps in (1.7) is a direct consequence 



of Proposition 4.7, (see also Remark 5.3). □ 



Remark 5.2. Note that M{Agx) generally contains high dimensional faces and therefore the chain 
expansion of Xrei '^x[gi,wj] involves monomials in eg and Wr- To simplify this expansion one 
may observe that M{Aqx) has a homotopy type of finitely many points or is empty. Specifically, 



from (4.5) we have 

XreKA,Aax)=x(A)-#{AnaX}. 
The random variable #{A n dX} (counting points in Aqx) can be expressed as follows: 

g 

XreKA,Aax) =X(A)- J^li^li}. (5.5) 

i=l 

where {1,. . . ,q} label points of dX and {w{i}}j=i,...,5 are the indicator functions of points in A^x- 
Consequently, we may derive expressions for powers xtei polynomials in ]R[e/, u){jj,]. These 
expansions of Xrei involve products of Cg and w^ij only, which may provide a different way to 
express F{X C\A\). 

Remark 5.3. In order to be more explicit about how the computation of ps,r simplifies in the 
case the nerve AA(A) equals the Vietoris-Rips complex Tl{A): let us suppose Aj are e-radius closed 
balls in X with random centers G X. In 7^(A) any simplex indexed by I = {ii,i2, ■ ■ ■ ,ik} is 
determined by its edges, and an edge {i,j} in 7^ (A) occurs if and only if — < 2e (where 
I • — • I is a short notation for the distance •) on X). For instance, we have 

PI = P(Aii n A,2 n . . . n Ai, / 0) = p(|ei, -ii,\<2e\ ys,t s^t). 



Enumerate points in dX as follows {xi,X2, ■ ■ ■ , xm}, M = j^{dX}. Now, ps^r given in (5.1 ) is just 
a volume of the set 

-4s,r = {{£.1, . . . ,^„) G X" I V/esVs,te/, |6 - (.t\ < 2e, V/er 3i<s<m ^if^i |6 - Xs\ < e}. 

which in the case F = d^i d^2 ■ ■ ■ d£,n (i-e. ^i's are independent) can be computed via ordinary 
calculus techniques or estimated numerically. These formulas further simplify, if dX = 0, but we 



do not attempt these computations here, except in the test case considered in Section |7.1 

6. Upper bounds for coverage probability 
In this section we use the method of finite differences, c.f. pLj, to estimate the complete coverage 



probability from above in terms of the mean of the Euler characteristic and prove Theorem 1.2 
Let {Aj}, i = 1, . . . ,n be a finite good covering of X, consider the following shifted version of the 
relative Euler characteristic Xre/(A, Agx) of {M{A),M{Aqx))- 

Xo =Xrei{A,AQx) -m, 



where m = Xreli^^ dX). From (2.9) we obtain 

Xrel{A,Adx)=Po-Pl, (6.1) 

where /S* = j3*{A, Aqx) stand for the random relative Betti numbers. Recall that {ej,wj}, 
f{n) stand for the indicator functions of faces in {M {A) , M (Agx)) ■ 
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We will consider a filtration by random vectors Vi denoting (e7(j), where J{i) E f{n) 
are subsets of {1, ... Note that Vi reveals subcomplexes in spanned by vertices 1 through 
i. By analogy to the setting of Erdos-Renyi model [Tj, we set up a vertex exposure martingale, 
associated with xo and {Vi} as follows: 

Yo = Ho = E{xo), Y^ = E{xo\V), i = l,...,n, (6.2) 

Clearly, = Xo and the sequence {Yj} is an instance of Doob's martingale |Tj. Recall the following 
variant of the Azuma-Hoeffding inequality [Il[3], for {Yj}: 

2 

- Yo < -a) < exp( ^ ~" (6.3) 



where a > 0, and Cj is a difference estimate 

|Yi -Yi_i| < Ci. (6.4) 

Exposing a vertex (or a face containing it) changes ySo by at most 1 and by at most Pi{X, dX) = 
1 — Xrei{X, dX) thus we obtain 

|Y.-Yi„i| <2 + \Xrel{X,dX)\. 

Let a = /iQ, then 

P(Xo = 0) = P(xo < /io - a) = P(Y„ - Yo < -a). 



Using the above estimates for and (6.3) yields 



P(X C |A|) =P(xo = 0) <exp' 



2n{\xrei{X,dX)\ + 2Y)' 



which completes the proof of Theorem 1.2 



7. Distribution of x on S"" and Steven's results. 



In [3T] Stevens proved Formula (1.1) for the complete coverage probability by independent identi- 
cally distributed n random a-length arcs on the unit circumference circle X = . Naturally, it is 
a case of a random covering A = {Aj}, i = 1, . . . , n on X = we consider in this paper. Work in 
|31j provides more information, i.e. exact distribution of number of gaps G = G(A) = 0oiX — |A|) 
in X — |A|. It is given as follows: let k = min([^J,n) then |3H p. 319] 



P(G(A)=i)= Q^(-l)'-^(^_^')(l-ia)'^-\ for 0<j<k 



P(G(A) = j) = 0, for j > k 



(7.1) 



where the coverage probability equals P(G = 0) (i.e. no gaps). By Proposition 4.7 the covering 
A is good provided a < |. The random Euler characteristic of the nerve AA(A) (thanks to (2.9)) 
satisfies 

x(A)=;9o(A)-;9i(A) 



where (thanks to the Nerve Lemma 4.1) the random Betti number l3o measures the number of 
connected components of A and /9i(A) a number of cycles in A, where the later in the case of 
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cannot exceed 1. If y8i(A) = we have = /0o(A). Because any realization of |A| is just a union 
of arcs it follows that 

X(A) = G(A), (7.2) 



Thus, (7.1) is also the distribution of x(A). In the following subsection we compare formulas in 
(7.1) and those obtained from Theorem |1.1| in the case of three arcs (n = 3). 

7.1. Verification of the Steven's formula in the case of three arcs. Suppose we are given 
three arcs {A, B, C} on the unit circle of length a and suppose that arcs' centers {X, Y, Z} 
are distributed uniformly and independently on S^. In order to simplify calculations we will just 
compare the moments of x(A) obtained from Equations (1.12) and (1.17) to the moments of G := 
G(A) directly computed from the distribution in (7.1). In the relevant case of three arcs of length 
a < 1, possible values of k = min([^J,n) are 2 and 3. Formula (7.1) yields the following 

'3^ 



Po = P(G = 0) = 1 



{i-ay + 



[1 - 2ay 



-(1 - 3af 



1? or 







Pi 



P2 



1 



a 



(3a 
3 

3+ 18a 
3 

3(1 - 2a)^ or 3a(2 - 5a) , 



21a^ 



2 
1 

or 



l-2a) +(l-3a)^ 



6a^ 



i) = (2) ((1 - 2a)^ |-(l-3a)^| ) 



P3 = P(G = 3) = +(1 - 3a) 



where the "unboxed" expressions are valid for k > 2 and "boxed" expressions for k 
we obtain direct formulas for the moments 

2 



3. As a result 



E{G) 



: Pi + 2p2 + 3p3 
: Pi + Ap2 + 9p3 



E{d^) =pi + 8p2 + 27p3 



3 — 6a + 3a 
9 - 30a + 27a^ 

_ f21 - 78a + 75a2, 
1 27- 114a + 129a2, 



for 
for 



k 
k 



2, 
3. 



(7.3) 



We aim to verify the above by applying the technology of Theorem |1.1[ As a first step we compute 
the discrete probability measure defining the random nerve complex A/'(A). Let us denote by x, 
y, z the characteristic functions of vertices (i.e. ejjj's), then o, b, c the indicator functions of the 
edges (i.e. e^j jj's) and F the indicator function of the face (i.e. Cji 2,3})- With this notation, the 
random Euler characteristic is represented in Mx[x, y, z,a,b,c, F] as 



X=x+y+z- 
= X + y + z — 

From the underlying assumptions 

P(A = 0) = P(B = 0) = P(C = 0) = 



a-b-c+F 

axy — byz — cxz + Fabcxyz. 



thus P(x = 1) 



1) 
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From (1.17) we obtain 

x'' = {x + y + z) + {-2 + 2^){xy + xz + yz) + {?, - 2, ■ 2^ + 3^)(xyz) 

+ (1 - 2^)(axy + byz + cxz) + (-1 + 2^"+^ - 3'')(axyz + hxyz + cxyz) 
+ (1 - 2^+^ + ?,'^){ahxyz + hcxyz + acxyz) + (-3 + 3 • 2'' - 'i^){abcxyz) 
+ Fabcxyz, 

and after substituting x = y = z = 1 

X = 3 — a — b — c + abcF, 

= 9 - 5{a + b + c) + 2{ab + ac + be) + a6cF, 
X^ = 27- 19(a + 6 + c) + 12(a6 + ac + be) - 6abe + abeF. 



(7.4) 



(7.5) 



Next, we need distributions of the Bernoulli monomials: a, b, c, ab, ae, be, abe, abeF which represent 
all possible realizations s G of AA(A). Recall that the arcs centers are uniformly distributed on 



V 



V = U + {l-a.) 



1-a 




V = U -a 



V = U -{\-q) 



U 



Figure [7?T| Regions on the unit square representing P(afoc) and P(a6c_F'). 



and independent, we denote the distance on 5^ by | • |. Due to the symmetry it suffices to 
compute only probabilities for a = axy, ab = abxyz, abe = abexyz, and F = Fabexyz. 

¥{axy) = F{{\X -Y\<a}) = [[ 1{{\X -Y\< a}) dX dY dZ 

dUdVdW = 2a. 

'{\U\<a} 

where in the third identity we applied U = X — Y , V = X , W = Z . Analogously, we have 
W{abxyz) = ¥{{\X -Y\<a,\X - Z\< a}) 

m\U\ <a,\V\ < a}) dU dV dW = Aa'^ , 
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F{abcxyz) = F{{\X -Y\<a,\X - Z\< a, \Z-Y\< a}) 

H{\U\ < a,\V\ <a,\U-V\< a}) dU dV dW 



_ JSa^, " < i 

\3a2 + (3a- 1)2^ «>i 

for U = X — Y,V = X — Z,W = Z. The last computatfon is represented on Figure |7.1| which 
pictures one of the faces of the torus (5*^)^ as a unit square in U, V coordinates. F{abcxyz) is the 
sum of areas outside the hnes V = U ziz {1 — a) and inside the hues V = U ziz a, i.e. the sum of 
the checkered pattern region and the cross pattern region. These regions have a direct geometric 
interpretation the checkered region corresponds to configurations of arcs which have non-empty 
triple intersection and the cross region of area (3a — 1)^ corresponds to configurations of arcs which 
have only pairwise intersections (these are the configurations which actually give us coverage of 
S^). Note that the checkered region disappears for a < ^ which is consistent with the fact that in 
this case it is impossible to cover the unit circle with three arcs. 

Computing ¥{F) directly is a bit more involved, it requires finding the volume of {| "^^^ — Z\ < 
% + 1 I } in the torus. We omit this calculation and use the observation of the previous paragraph 



to conclude that ¥{F) equals the area of the cross pattern region of Figure 7.1 and thus 

F{xyzabcF) = 3a^. 

Collecting the above information yields the probability distribution of the random nerve of {A, B, C} 
as follows 



• 

• 

• 


F{x) = P(y) = 
¥{z) = 1 


» 


r{axy) = fibyz) = 
r{cxz) = 2a 




¥{abxyz) 

— ¥{acxyz) = 

¥{cbxyz) = 4q2 


4 


F{abcxyz) — 3a^ 
for a < |, or = 
12a^ - 6a + 12 if 




r{Fabcxyz) = 3a^ 





As the last step we compute moments of x, the set of possible realizations of x is given as {0, 1, 2, 3}, 
three first moments suffice to recover the probability distribution of X- After substitution we obtain 

E{x) = 3 - 3F{axy) + F{xyzabcF) = 3 - 6a + 
E{x^) = 9 - 5P(axy) + 2¥{abxyz) + F{xyzabcF) = 9 - 30a + 27a^ 
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E{x'^) = 27 - 19F{axy) + 12F{abxyz) - 6F{abcxyz) + F{xyzabcF) 

_ J21-78a + 75a2, a < 5, 
~ \27-114a + 129a2, a > 



These results are in full agreement with those obtained in (7.3) 
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